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Abstract 



We study the effect of soft gluon emission in the hadroproduction of gluino-gluino and squark- 
antisquark pairs at the next-to-leading logarithmic (NLL) accuracy within the framework of the 
minimal supersymmetric model. We present the calculation of the one-loop soft anomalous 
dimension matrices controlling the colour evolution of the underlying hard-scattering processes. 
The numerical results for resummed cross sections for proton-proton collisions at the Large 
Hadron Collider are discussed in detail. 
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1 Introduction 



Supersymmetry (SUSY) is one of the most promising candidates for the theory of physics beyond 
the Standard Model (SM). In the coming years, experiments at the Large Hadron CoUider (LHC) 
will undertake searches for new physical phenomena. A large part of this effort will be devoted 
to looking for signals of SUSY. 

One of the most studied SUSY models is the Minimal Supersymmetric Standard Model 
(MSSM) [1], characterized by the minimal content of supersymmetric particles and i2-parity 
conservation. Within the MSSM, the dominant production processes of sparticles at the LHC 
are those involving pairs of coloured particles, i.e. squarks and gluinos, in the final state [2]. 
The exact discovery reach of the LHC is model-dependent but it is expected that the discovery 
of squarks and gluinos should be possible for masses of up to around 2 TeV [3]. Since the 
hadroproduction cross sections for strongly-interacting sparticles depend only on the masses 
of squarks and gluinos [4, 5, 6, 7], measurements of total cross sections for coloured sparticle 
production may be used to determine values of the fundamental MSSM parameters, e.g. the 
masses of sparticles [4, 6, 7, 8], or to draw exclusion limits for the mass parameters [9, 10]. The 
precision of the mass determination (exclusion limit) will crucially depend on the accuracy of 
the corresponding theoretical predictions. With the large production rates expected at the LHC 
it is thus of utter importance to study the total cross sections for the hadroproduction of squarks 
and gluinos with the highest possible theoretical accuracy. 

The leading-order (LO) total cross sections of 0{a^) were calculated long time ago [4, 5]. The 
corresponding next-to-leading order (NLO) SUSY-QCD corrections are known for all hadropro- 
duction processes of pairs of squarks and gluinos [6, 7]. The NLO electroweak (EW) correc- 
tions [11] for the processes involving squarks in the final state are also known, as well as the LO 
EW O(a^) total cross sections and the 0{aas) LO EW-QCD interference predictions [11, 12]. 

The NLO SUSY-QCD corrections have been found to be positive and large. Among the pair- 
production processes of coloured sparticles at the LHC, the gluino-pair (gg) production receives 
the largest NLO SUSY-QCD correction that may reach 100% for gluino mass rrig = 1 TeV and 
1.2 TeV < niq < 2 TeV [7]. The corrections to the squark-antisquark (qq) total cross section 
can be also sizable, of order of 30% for the squark mass niq = 1 TeV, and are the second 
largest in a certain range of mass parameters. The occurrence of large corrections indicates 
that computation of higher order SUSY-QCD corrections is necessary in order to achieve precise 
theoretical predictions. 

A large part of the NLO SUSY-QCD corrections to the total cross section for gg and qq 
processes comes from production close to threshold [7]. The threshold region is reached when 
the square of the partonic ccntcr-of-mass (c.o.m.) energy, s, approaches 4m^, where m is the 
average particle mass in the produced pair. The velocity of the produced heavy particles in 
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the partonic c.o.m. system (3 = \/l — Arri^/s is then small, (3 <^ 1. In this region two types of 
corrections dominate: Coulomb corrections due to exchange of gluons between slowly moving 
massive particles and soft gluon corrections due to emission of low energy gluons off the coloured 
initial and final states. The soft gluon corrections are enhanced by powers of large logarithms 
of (3, i.e. at the NLO one finds terms in the relative corrections proportional to ag log^(/3^) that 
become sizeable when 0^ ~ exp(— l/y^). At the n-th order of the perturbative expansion in 
the strong coupling CKg the total cross sections receive corrections proportional to a" log'^(/3^) 
where k = 2n, . . . , 0. Sufficiently close to the partonic threshold fixed-order expressions for 
the cross sections are bound to fail. However, the logarithmic contributions can be taken into 
account to all orders in ag by means of threshold resummation. Resummed predictions are 
particularly important for processes with large masses in the final states since then the bulk of 
production comes from the threshold region. This is exactly the case for production of sparticles 
which are expected to be heavier than the SM particles. Additionally, if partonic subprocesses 
involve gluons in the initial state, the soft-gluon effects, and thus the impact of resummation, 
are expected to be significant due to the high colour charge of the gluons. 

Crucially, calculation of the soft gluon corrections provides a reliable estimate of unknown 
higher order terms beyond the NLO, what results in reduction of the theoretical uncertainty 
due to scale variation. In a recent letter [13] we have presented results for threshold-resummed 
cross sections at the next-to-leading logarithm (NLL) accuracy for hadroproduction of gg and 
qq pairs at the LHC. A dominant part of the next-to-next-to- leading (NNLO) correction for the 
qq production consisting of terms coming from the expansion of the resummed exponent at the 
next-to-next-to-leading logarithm (NNLL) level. Coulomb terms and the universal scale terms, 
was later calculated in [14]. Moreover, threshold resummation for single colour-octet scalar at 
the LHC was also studied [15]. Further work on resummation for production of coloured sparticle 
pairs is to be found in [16, 17]. 

The hadronic production processes of gg and qq pairs are scattering processes with a non- 
trivial colour flow structure. At the NLL level resummation requires including contributions from 
soft gluons emitted at wide angles. Such emission is sensitive to the colour flow of the underlying 
hard scattering and the evolution of the colour exchange is governed by the soft anomalous 
dimension matrix [18, 19, 20, 21, 22]. The one-loop soft anomalous dimension matrices were 
first calculated for heavy-quark and dijet production [18, 19, 21, 23]. In [13] we have presented 
the explicit form of the one-loop soft anomalous dimension matrices for partonic subprocesses 
contributing to gg hadroproduction. The general results for any 2 ^ n QCD process with 
massless particles in the final state were derived at one- [22], and two-loops [24]. The two- 
loop anomalous dimension for the pair production of heavy quarks was also determined in the 
threshold limit [25] , and employed later in [14] . Recently, the structure of the massive two-loop 
matrix for any 2 ^ n process has been studied in [26] . 
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In this paper we discuss in detail the derivation of the analytical results presented in [13] 
and carry out a thorough study of the numerical results as well as perform resummation of the 
leading Coulomb corrections for the qq and gg production processes. 



2 Leading order results 

The hadronic cross section for the process /11/12 — > kl reads 

(^hihi^kliSjirri^}) =J2 dxi dx2 fi/h-^{xi, jip) fj/h2{x2, fJ-p) o-ij-^fe;(s, {m^}, , (1) 

where S {§) is the square of the hadronic (partonic) center-of-mass energy and {m?} stands for 
all masses entering the calculations. The parton distribution functions (pdfs) fi/h are taken at 
the factorisation scale /xp. We set equal to the renormalisation scale fiR in our calculations. 

The expressions for the LO partonic cross sections aij^kl{s, {m^}, /xl^, //^) for all squark and 
gluino hadroproduction processes can be found in [7] . Here we present the contributions to the 
LO cross sections for the qq and gg production coming from different colour channels. For each 
partonic process we define the corresponding colour basis in the s-channel. 

For the qq production we consider the processes 

qi{Pi,Oii) q{pj,aj) q{Pk,Oik) q{PhOii) (2) 

and 

9{Pi,ai)9{Pj,aj) ^ q{Pk,ak) q{Pl,Oil) , (3) 

where p are particle four-momenta and a and a are colour indices in the fundamental and adjoint 
representation of SU(3), correspondingly. In the quark-channel (2) we have only two possible 
colour exchanges: the singlet and the octet, {1,8}, and the basis consists of two colour tensors 

In the gluon channel, the basis is built out of three tensors corresponding to {1, 8s, 8a} repre- 
sentations 



c£'5 = Tt^. d 



b rihai a. 



8s OilOLk 

41 = Kaj''''''' ■ (5) 



where T'^ matrices are the SU(3) generators. 
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At the leading order, two partonic channels contribute to the gg production: 

q{Pi, "i) Q{Pj,o:j) -> 9{Pk, ak) g{Pl,ai) , 



(6) 



and 

9{Pi, ai) 9{Pj,aj) g{Pk,ak) 9{Pl, ai) ■ (7) 

For the process (6) the colour basis is the same as for the qq production in the gluon-channel (3) 
and is given by (5) after interchanging the indices (ij) <r^ (kl). For the gg channel there are eight 
independent colour tensors. Following [21] we choose an orthogonal basis, {cf '^}, / = 1, 2, . . . , 8, 
consisting of five tensors cf'^, cf'^, cf'^, cf'^ and cf'^ corresponding to the 1, 8s, 8a, 10 © 10 
and 27 representations in the s-channel, and three additional tensors, Cg'^, cf'^, and cf'^. The 
base tensors are 



^3 



^,9 



^6 



4'^ 



5 ' 

3 



2 8 



^Uiajb^bakai 



^ ^ j!aiajb^bai;ai _|_ ^aiUjb j^bakai"^ 
^ ^jaiQjb^bakai _ ^aiajb jbakai^ 
— ^(J^"-i"-kb j^bajai _|_ j^aia^b ^baj ai'^ 



(8) 



In the set of basis defined above, we obtain the following colour-channel contributions to the 
total cross section for qiqj — ^ qq: 
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gg-*qq,8s 
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and for gg — > ^5: 
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(20) 



(21) 



/^ .§(l-/3,-) + 2m2_ ' 
= l,s(l + /3,-) + 2m2__ 

where as is the SUSY Yukawa coupling. 



Lo = loe 



's(l-/?3)-2m2 
s(l + p-g) - 2m? 



3 Threshold resummation — general framework 

The resummation for 2^2 processes with all four external legs carrying colour was studied 
extensively in the literature. The resummed cross section for the heavy-quark production was 



5 



constructed in [19, 23], and for the dijet (multiple jet) production in [20, 21, 22]. Here we briefly 
review the derivation of the resummed cross sections for the production of two coloured final 
state particles of equal mass m, in the form presented in [23]. In our calculations me make use 
of the framework of [19, 20, 21]. 

Using the hadronic threshold variable p = 4m^/S we rewrite the cross section (1) 

(Thih2-^klip, {'m'^})=^Yl / dxidx2 dpS(p — ) fi/h^ {xi , p)fj/h2 ix2, p) (Tij-^kl,iip, {m'^}, 1^'^), 

J J \ X1X2J 

(22) 

where the index / sums over all possible colour states of the hard scattering. 

At higher orders in perturbation theory, the partonic cross section a contains terms of general 
structure a" log™ m < 2n, with f3 = \/l — p. These terms arc singular in the threshold limit 
p ^ 1. The singularities can be systematically treated by taking Mellin moments of the cross 
section 

^hih2-*kl{N,{m'^}) = I dp p^~^ ahihi-^klipdrri^}) (23) 
Jo 

i,j 

The moments of the parton distributions fi/h{xi,p'^) are defined in the standard way, ^ 

Lh{N,p^)= ['dxx''-^h/h{x,p^), (24) 



i/hV^,P- ) = 'J'XX Ji/h 

and the moments of the partonic cross section ab — > kl are given by 

aij^kliN,{m'^},p'^) = I dpp^^^ aij^kl{p,{^^}^P^) ■ (25) 

Taking the Mellin moments transforms the logarithms in 0^ into the logarithms of the Mellin 
variable N which are then resummed to all orders in CKg. 

Following [27] we define the differential distribution 

^^^g^^' (p,{mn) = (26) 

= XlZl / dXldX2 fi/h^{xi,p) f)lh^{x2,p)5{^- XiX2)(Tij^kl,l{pl^^{^'^)^l^'^) 
1,3 I 

rC'+ioo 

= EE^^J-fc^/W^'W}'/^') / a!iV ^^/^//.l(iV,M') (iV,/x2). 

i,j I JC-ioo 

Since soft radiation carries colour charge, it can change the colour state of the underlying hard 
scattering for hadronic processes with two or more coloured partons in the final state. This 



^Note that from now on we will use the tilde sign to mark symbols for iV-space quantities. 
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has to be taken into account while writing the form of the cross section with long-distance and 
short-distance effects factorised, and leads to [20, 23] 

^''r'\ pdm'}) = EE h*j^kl,i{p/^,{m'},fi')hij^kl,j{p/^dm'},fi') (27) 

i,j l,J 

Y fC+ioo 

X 



J C—ioo 



2TTi Jc 

with = 4m?. The function hij^klj ( j) is a colour-dependent hard-scattering ampli- 

tude (conjugate of) absorbing the far, i.e. of the order of the scale of the process Q, off-shell 
effects. All the logarithmic dependence on N, originating from soft and collinear radiation, is 
contained in the function oJ^j'^^j^i Q,fi)- 

In the approach of [19, 20, 21, 28], resummation follows from refactorisation of partonic cross 
sections. In the case of threshold resummation, the cross sections are factorised w.r.t. fixed 
fractions of energy as opposed to fractions of momenta in the standard collinear factorisation. 
Using the refactorised form of the cross section for the production of two massive coloured 
particles in the final state [19], we have, up to corrections of 0{1/N), 

'^''-''-{p,{m'}) = Y.T.h:j^klAP/^,{m'},n')h,j^M,j{p/^dm'},n') (28) 



2TTi Jc 



C+ioo 
C—ioo 



fi/i{N, fj,^)fj/j{N, n^) 

Following [19, 20] the soft eikonal function S represents coupling of the soft gluons to the initial 
and final state particles. Consequently, the soft function carries dependence on the possible 
colour exchanges and we sum over all possible colour structures /, J at hard vertices. As an 
object of purely eikonal character [19, 28, 29] the soft function can depend on the scales only 
through their ratio. The parton-in-parton distributions V'i/j are defined at fixed fraction of energy 
of parton i in the partonic center-of-mass frame as opposed to the light-cone distributions /j/j 
which are defined at fixed momentum fraction. The parton distribution functions and the soft 
function can be defined explicitly in terms of operator matrix elements [20] . 

As a consequence of refactorisation, the soft function Sijjj and the distributions functions 
V'i/j) fi/i obey the corresponding rcnormalisation group equations (RGBs) [19, 20, 21, 28]. So- 
lutions of these RGEs give functions which rcsum the large logarithms in question. In fact for 
the 2^2 production process involving coloured massive particles only the solution of the RGE 
for the soft function 

P'-qJ^ + ^(S-)^ j Sij^kl,IJ = -'^ij^kl,IK^ij-*kl,KJ - Sij^kl,IL^ij-^kl,LJ (29) 
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is needed. As shown in [19], the resummed initial-state jet factors for the ratio of ^/ij/j to /j/j 
functions in Eq. (28) can be obtained directly from the resummed cross section for production 
of a colour singlet state through the Drell-Yan mechanism, for which results are known [28, 30]. 
We write this ratio in the form 



r n —1/2 

Ri{a,{^?))^i{N,Q\^?) U,-{Q/{Nti),as{fi^)) , (30) 



where the function Ri{as{fJ?)) is an A^-indepcndcnt and infrared-safe function of the coupling. 
The soft cikonal function Uf^ describes soft gluon emission and exchange by the annihilating 
initial state partons in the Drcll-Yan process. The radiative factor Aj represents both the soft 
and coUincar radiation from an incoming parton. Inserting the expression for the ratio of to 
/ functions, Eq. (30), into Eq. (28) leads to 

d ^^^^^ 

= Y.Y.hh^kl,iip/^,{rn'},fi')h,,^kl,j{p/^,{rn'},^i') 

"I rC+ioD 

X TT- / dNr^f,/h^{N,fi^)f]/f,^{N,fi^)Ri{a,ifi''))Rj{a,{^i^)) 

ZTTt J C—ioo 

X Ai{N,Q\i?)Aj{N,Q\i?)~Sij^kijj{Q/{Nij,),a^{ix^)) , (31) 

where we introduce Sij^kljJ = Sij^kl,lj/U{i. Consequently, the soft anomalous dimension 
matrix T^j^j^i corresponding to the function Sij^kljJ is given by 

^ij-*kl,Ij{ois) = ^ij-*kl,ij{ois) - Slj'^iiiois) , (32) 

where Fj^ is the anomalous dimension associated with the Drell-Yan soft function U{i. The 
results for T^j can be found in the literature [21, 28]. 

In general, for a given colour basis, the soft anomalous dimension matrix Fjj is not diago- 
nal, leading to resummed expressions in terms of path-ordered exponentials [20]. Through the 
diagonalisation of the soft anomalous dimension matrix a simpler form for the resummed cross 
sections, involving a sum of exponentials, can be obtained. However, the diagonalisation proce- 
dure can be avoided if from the beginning the calculations are performed in the colour basis in 
which the the soft anomalous dimension matrix is diagonal, i.e. 

Tij^klias) = diag(. . . , Xij^kljias), • • •) • (33) 
The soft function then reads, up to NLL, 

,2^^ _ q(o) 

kl, 



X exp 



/ U*j-.jfc/,/(Q!s(9^)) + Kj^kl,j{as{q'^)) 



, (34) 
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with the 0-th order term in the perturbative expansion of Sij^kijj{l,as{Q'^ /N"^)) straightfor- 
wardly related to the colour structure, i.e. 

'Stlkljj = Tr(4c7) , (35) 

where {cj} is the corresponding colour basis for the process ij kl, see Section 2. If additionally 
the colour basis is orthogonal then the soft function matrix becomes diagonal 

rQ/N dq 



/ -^2Re(A,,-^«,/(as(g'))) 



In this case inserting the solution (36) into Eq. (31) leads to the following NLL expression 

1 /-C+ioo 



(36) 



where we identify 



I hij-,ki,i I ^ RiRj sfjL,ki,ii = (^ijLki,Aj-^ki,i ■ (38) 

The functions Cij^klj a^e of perturbative nature and contain information about higher-order 
corrections which are non- logarithmic in TV. All the information on the soft non-coUinear log- 

(s) 

arithmic corrections is included in the radiative factor A-i^j^i j. After performing integration 
over ^ we obtain 

i,j I JC-ioo 

X A,{N,Q\^^^)Aj{N,Q\^?)A^_^^^J{N,Q\^''). (39) 

The functions Cij^kij a-re related to the Cij^kij functions in Eq. (37) and have a perturbative 
expansion of the form Cij^kij = 1 + J2n=i ^s^ij^kl i- general, the values of the coefficients 

-(n) 

^ij-^kl I obtained by comparing the resummed cross section expanded in with the ex- 
pression for the full higher-order cross section in N space. The expression for the resummed 
hadronic cross section in the Mellin-moment space can be easily derived from Eq. (39) and reads 

^Ti^kli^^ {"^'}) = E E h/H, (iV + 1, /x^) (iV + 1, ag-i,, ,(iV, {m^}, m^) 

X Cij^kiA^. {m'}, /x^) Ai{N + 1, Q\ fi^) Aj{N + 1, Q^, /x^) aJ^^, ^(A^ + 1, Q^, fi^). 

(40) 

(s) 

The expressions for the radiative factors Aj, A^^i^^^ ^ arc presented below. Since after expansion 
of the exponentials the non-trivial terms contained in C^j^kij generate contributions of the 
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NNLL and higher orders, we keep Cij^^ij = 1 for the rest of the calculations. Additionally, 
in Appendix A we list the results for the Mellin moments of the colour-channel contributions 
^ij^-^kl I *° leading-order partonic cross sections. 

3.1 Soft radiation factors 

The expressions for the radiative factors in the MS factorisation scheme read, up to the NLL 
level [19, 20, 23], 

logA,(Ar,Q2^/x2) = dz J ^ ^A,{a,{q^)), 

As already noted, the radiative factor Aj represents the soft and coUinear gluon radiation from 

(s) 

the incoming partons, whereas the function j takes into account soft and large-angle 

gluon radiation. The coefficient Ai is a power series in the coupling constant Og, Aj(as) = 
^Ap-^ + {^f A/^) + The universal leading logarithm (LL) and NLL coefficients \ Af^ 
are well known [31, 32] and given by 

.f).ic.((|4)c.-^„,). ,42, 

with Cg = Ca = 3 for radiation off gluon lines and Cq = C-p =4/3 for radiation off quark lines. 
In the case the soft anomalous dimension matrix Tij^^i is diagonal in the orthogonal colour basis, 
the relation between the -D^jl,;,; ^ coefficients and the Tij^kl eigenvalues, Xij^kl,i-> reads [19, 20, 
21] 

4-«,/ = 2Re(A,,^fe^,/(as)). (43) 



as 



TT 



The customary NLL expansions of the radiative factors (41) are presented in Appendix B. 



4 Soft anomalous dimension matrices 

Due to the same colour structure the soft anomalous dimension matrices for the qq production 
and for the heavy qiiark production are the same. The results for the heavy quark production 
are available in the literature [19]. For the gg production, however, a separate calculation is 
needed. 

In order to obtain resummed cross sections up to the NLL accuracy we need the one-loop 
0{as) result for the soft anomalous dimension matrix, r> 1^^.^ jj . The T 

ij^ki.ij matrices are 

given in terms of matrices of renormalisation constants Z^j^i^ijj for the soft function [18]-[21], 
[29] . In the MS scheme in the d = 4 — e dimensions we have 

^\]lki,iji9s) = -|^Res,^o^^fc;,,j(5., 6) . (44) 



10 



(a) 




ib) 



Figure 1: Vertex (a) and self-energy (b) one-loop contributions to the soft function Sij^kljJ 
with i,j massless and k,l massive particles. 

The calculation of T^^jl^^; jj reduces then to evaluating the UV-divergent part of the one-loop 
correction to the Sij^MjJ function. 

In Fig. 1, following [19], we show the set of eikonal one- loop diagrams contributing to the soft 
function Sij^kijj- All one-loop integrals needed for the calculation of the anomalous dimension 
matrices for qq gg and gg — > gg are known. In our calculation we use the results for the eikonal 
one- loop integrals from Ref. [18, 19] which have been obtained in the axial gauge. Although 
the integrals are the same as in the heavy-quark production case, their contributions to the 
soft anomalous dimension matrices come weighted with different colour factors. We compute 
these colour factors in the set of basis presented in Section 2. The calculations are performed 
in two independent ways: using the FeynCalc package [33] for Mathematica, and using the 
Group Theory (Colour) Factors of Feynman diagrams package [34] for FORM [35]. The resulting 
expressions for the soft anomalous dimension matrices for the gg production are presented below. 
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4.1 Gluino-pair production 



We introduce the notation 



T = loi 

U = lo; 
S = - 



m 



t 



ITT 



V rn^s 

2 



m 



u 



2 

1 — ZTT 



(45) 



where the Mandelstam variables are given by 

and = i(l - 2mVs) (log + i7r). We also define A = f + U, n = f-U. 

In the basis (5) we obtain the one-loop soft anomalous dimension matrix of the form 



Qq^gg 



55(«s) 



L V 



65 



\ 



35+|A 



20 
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-ZTT I 



(46) 



The one-loop soft anomalous dimension matrix for the gg channel, calculated in the basis (8) 
has the block form 



f« -( 
gg^gg^ 

where the five-dimensional matrix f 5 reads 
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and the three-dimensional matrix T3 is diagonal, 

fa = diag(3(5 + C7), 3(-S + r), 3(f + £/) ). 



(48) 



(49) 



Although both Sij^kl,l.J a-nd Uf^ arc gaugc-dcpcndcnt functions, results for the one-loop 
-^ij^fc/ // matrices presented here arc gauge-invariant. The gauge dependence in the sum of ver- 
tex and self-energy corrections cancels against the gauge dependence of the Drcll-Yan anomalous 
dimension [19, 36]. The U^i soft function, in turn, carries the same gauge dependence as the 
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ratio of the incoming jet functions {i^i/il fi/if [21, 29]. In practice, we fix the gauge AP = Q to 
calculate T^j^j^.^ jj and use r^l^(as) = Cia^j-K in Eq. (32). 

At the production threshold Pg ^ the soft anomalous dimension matrices Vqq^gg and 
Tqq^gg approach the diagonal form. The off-diagonal terms, proportional to fl, vanish like Pg 
for Pg and thus may be neglected. Using Eq. (43) we obtain 

{DS^gg,l} = {0,-3,-3} (50) 

{^£1^5, /} = {0, -3, -3, -6, -8; -3, -3, -6} , (51) 

where the index / indicates a colour channel, defined by the base tensor cf^ and cf^ as in (5) 
and (8), correspondingly. Note that the values of the D'^-'^^-coefficients are the negative values 
of the quadratic Casimir operators for the SU(3) representations for the outgoing state. This 
agrees with the physical picture of the soft gluon radiation from the total colour charge of the 
heavy-particle pair produced at threshold [23]. 



4.2 Squark-antisqueirk production 

For completeness, we also list here the NLL coefficients D^^^ , which we need for our numerical 
calculations. They have been first obtained in the calculation of resummed cross sections for 
the heavy quark production and read [19, 23] 

{^SUv> = {0>-3} (52) 
{^SU7} = {0'-3>-3} (53) 

with the index / indicating the corresponding tensors in the s-channel colour basis Cj'^ and c|'*, 
given in (4) and (5). 



4.3 Checks of analytical results 

The computational framework applied here has been tested by re-deriving the known results 
for the one-loop soft anomalous matrices for the pair-production of massive particles. We have 
reproduced the results for the heavy quark production [19] both in the qq and in the gg channel. 

Another test of our calculations has been based on a comparison between expansion of the 
resummed cross section, Eq. (40), and the Mellin moments of the NLO corrections taken in the 
threshold limit. The analytic form of the NLO corrections in the threshold limit in momentum 
space is known for all squark and gluino production processes [6]. More precisely, for each of the 
partonic processes, qiqj qq, gg qq, qq gg, and gg gg, we have extracted terms with 
log^ and log A from the 0{as) expansion of the corresponding resummed formula. In each 
case the result has been compared with the collection of terms logarithmic in A in the Mellin 
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transform of the 0{as) correction taken in the large N Hmit. The contributions enhanced by 
double logarithms, i.e. 0{as log^ N), depend only on the colour charges of the incoming partons 
and thus do not provide a check of the soft anomalous dimension matrices. A non-trivial cross- 
check is, however, provided by the contributions to the NLO correction with single logarithms, 
0{as log N), that are sensitive to the eigenvalues of the T matrices. Our result for the resummed 
cross section, Eq. (40), agree in this way with results of Ref. [6]. 



5 Resummation of leading Coulomb corrections 

Important higher order corrections to cross sections for production of coloured particles come 
from multiple exchanges of Coulomb gluons between the produced particles. This type of cor- 
rections should be then taken into account to all orders [17, 37, 38]. In the threshold limit, the 
Coulomb corrections are enhanced by the inverse powers of p. At n-th order in perturbation 
theory the leading corrections are of the form C^^^^j a" //?"■. These leading contributions can be 
summed to all orders using [37] 



;,iC) (0) 

where 



,A(^)(^..,^,,,), (54) 



A(^)(.) = , (55) 

exp(z) — 1 

and P = (3q for qq production, (3 = (3g for gg production. For the processes of interest, the k 
coefficients, calculated in the set of colour basis introduced in Section 2, are given by 

4 1 



qq^qq, I \ 3 ' g 

4 11 



^gg^qq,! \ 3' 6' 6 

- ( ^-^ A A 

'^qq^gg, ^ ~ \ ' 2 ' 2 

- ( --i _^ _^ n 1 _^ _^ 
'^gg-*ggJ ~ \ ' ~2 ' ~2 ~2 ' ~2 

The 0{as) Coulomb correction is a part of the full NLO result. Since wc arc interested 

in corrections above NLO, we subtract it from cjl-j^^^ j. In this way we define the correction 

(C) 

S^ij-^kl I leading (in terms of powers of 1//3) Coulomb contributions above NLO 



(56) 

(NLO) ^ ' 



'ij^kl,I — '^ij^kl,! '^ij^kl,! 
The corresponding Coulomb correction at the hadronic level then reads 

SafX^kliPA^'^}) = J2 dxidxidpSip —\ fi{xi,n) fj{x2,p)da^fl^ij{p,{m'^},n^) 

■ J. J J \ XXX2J 

(57) 
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6 Predictions for squarks and gluino production at the LHC 



We investigate in detail the effect of the soft gluon corrections on the cross sections for two 
sparticle production processes at the LHC, pp —>■ gg and pp —>■ qq, at y/S = 14TeV. The 
main results obtained in this section are the resummation-improved total cross sections. We 
also study the effect of the resummed leading Coulomb corrections to the total cross sections. 
Moreover, we present a detailed analysis of the soft gluon corrections in partonic channels of 
hadronic processes, including also the dependence of soft gluon effects on the colour structure 
of the hard matrix element. All numerical calculations were performed using two independent 
computer codes. 



6.1 Inversion and matching 

The resummation-improved cross sections are obtained through matching the NLL resummed 
expressions with the full NLO cross sections, 

-&SL(P,{-'}) = -fXijP^W'}) (58) 



■ ■ „ ^ „ CT 



1,J=<1,<1,9 



where is given through Eq. (40) together with Eq. (23), and cj^J'^^,; |(nlo) represents its 

perturbative expansion truncated at the order of ag associated with the NLO correction. 

The inverse Mellin transform (58) is evaluated numerically using a contour CT in the 
complex- TV space according to the "Minimal Prescription" method developed in Ref. [27] . More 
specifically, we use a contour parameterised by a parameter iV = Cq+x exp(ibz0), as described 
in [39, 40]. In order to be able to use available parameterisations of parton distribution functions 
in x-space we apply the method introduced in [39] . The NLO cross sections are evaluated using 
Prospino [41], the numerical package based on calculations employing the MS renormalisation 
and factorisation schemes. 



6.2 Numerical results 

In the phcnomcnological analysis we consider a wide range of gluino and squark masses. Left- and 
right-handed squarks of all flavours are assumed to be mass degenerate. For the gg production 
we vary the gluino mass, rrig, between 200 GeV and 2 TcV. Similarly, for the qq production 
we take ^ 200 GeV < rriq < 2 TeV. We present the results for a fixed ratio of gluino and squark 

^For the highest masses considered here, the experimental exploration at the LHC will require luminosities of 
0(100 fb-i). 
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masses, r = ^ and choose the fohowing values r = 0.5, 0.8 1.2 1.6, 2.0. The qq cross 

Tfiq 

section accounts for production of all qq flavour combinations apart from the ones with scalar 
top particles. 

For most of the phenomenological results, we use the CTEQ6M [42] parameterisation of 
parton distribution functions (pdfs). In addition, we give the total cross sections and study their 
scale dependence for the MSTW parameterisation [43] of the pdfs. Unless explicitly specified 
otherwise, the CTEQ6 pdfs are applied. In both sets of pdfs the usual assumption of five 
massless quark flavours active at large scales is made. Consequently, in the NLO and NLL 
calculations we use the two-loop MS QCD running coupling constant Og with nf = 5. We 
also show some results obtained at the LO accuracy using CTEQ6L1 parameterisation of the 
pdfs and the one-loop running coupling constant with 5 flavors. For all parameterisations we 
consistently use the corresponding default values of A^^) . The effects due to virtual top quarks 
and virtual sparticles in the running of as and in the evolution of pdfs are thus not included 
in our predictions. However, the value of the top mass mt = 175 GeV enters the matched NLL 
cross sections through the NLO corrections. 




0.2 0.5 1 1.5 2 0.2 0.5 1 1.5 2 

mg[TeV] m^[TeV] 

a) b) 

Figure 2: The NLL K -factor, JCnlL; for the gg (a) and the qq (h) total production cross section 
at the LHC as a function of gluino and squark mass, respectively; r = nig/mq. 

In Fig. 2 we present the relative enhancement of the NLO total cross sections due to soft 
gluon resummation, -R'nll = (j^™'^*'^^^/c7^^*~^. The NLL iC-factors are shown for gg and qq 
production cross sections at the LHC, in Fig. 2a and Fig. 2b, respectively. In the plots we set 
the scales = IJ-R = fJ-o, where fio = nig {fiQ = mg) for the gg production (the qq production). 
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The NLO CTEQ6M pdfs are used. -f^NLL grows with the final-state mass and depends on the 
mass ratio r in a moderate way. The relative correction -f^NLL — 1 reaches 16% (8%) for the 



gg production with r = 1.2 and rrig = 2 TeV (1 TeV), and 4% (2%) for the qq production with 
r = 2 and rriq = 2 TeV (1 TeV). The stronger effect found in the gg production follows from 
the dominance of the gg — > gg channel and hence larger colour factors. It comes from the fact 
that the size of the soft-collinear radiative factor Aj increases with higher colour charge of the 
incoming parton. Similarly, the size of the soft non-collinear gluon corrections increases with 
the total colour charge of the final state, which may be the highest in the gg gg case. 
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Figure 3: Scale dependence of the total gg (a) and qq (h) production cross section at the LHC 
(see the text for explanation). 

We also investigate the dependence of the matched NLL cross section on the values of 
factorisation and renormalisation scales, in comparison to the NLO cross section. To illustrate 
our results we choose fi = ^ip = fJ-R and r = 1.2. In Fig. 3a and Fig. 3b we plot the ratios 
^NLO(^ = ^^q)/^nlO(^ ^ and ^{match)(^ ^ ^;Uo)/fT('^^*=^) = fio), obtained by varying ^ 
between ^ = 1/2 and ^ = 2. Due to resummation the scale sensitivity of the gg production 
cross section reduces significantly, by a factor of ~ 3 (~ 2) at nig = 2 TeV {nig = 1 TeV). At 
nig > 1 TeV the theoretical error of the matched NLL gg cross section, defined by changing 
the scale fi = fip = fiR around /xq = nig by a factor of 2, is around 5%. In the case of the qq 
production the reduction of the scale dependence due to including soft gluon corrections in the 
theoretical predictions is moderate. 

The dependence on the pdf parameterisation of the NLL X-factors for gg and qq production 



17 




a) b) 

Figure 4: The relative NLL correction -JCnll for the gg (a) and the qq (h) total production cross 
section at the LHC as a function of gluino and squark mass, respectively; r = rrig/mq. The 
continuous lines correspond to the CTEQ6M pdfs and the dashed lines to the MSTW pdfs. 

at the LHC is presented in Fig. 4a and Fig. 4b, respectively. The MSTW pdfs, shown with the 
dashed Unes, lead to slightly larger NLL X-factors than the CTEQ6M pdfs. In this figure we 
choose r = rrig/mq = 1.2. The difference between the X-factors for the two parameterisations 
is moderate but it grows with increasing mass of the produced sparticles. It is expected since 
for larger masses of sparticles larger values of the factorisation scale and of the momentum 
fraction x are probed. For larger values of fip and x the pdfs are currently not so well constrained, 
see e.g. the discussion of the uncertainties of the pdfs given in [43] . 

Next, we analyse the combined dependence of the NLO and the matched NLL total cross 
sections on the factorisation and the renormalisation scale and on the type of the pdfs. The 
results are shown in Fig. 5 and Fig. 6 for the gg and qq production, respectively. Again, the 
renormalisation and factorisation scales are assumed to be equal, ^ip = h^R = ^ = '^A*0) where 
the scale /xq is set equal to the average mass of the produced particles. The parameter ^ is varied 
between 0.2 and 5. In Fig. 5a, Fig. 5b and Fig. 5c (Fig. 6a, Fig. 6b and Fig. 6c) the gluino 
(squark) mass, rrig (rriq), takes values of 0.5, 1 and 2 TeV, correspondingly, and r = mg/rriq = 1.2. 

The theoretical uncertainty of the pp — > gg cross section due to differences in the parame- 
terisations of the pdfs may be read out from Fig. 5. At the NLO the uncertainty is given by 
the difference between the dotted (CTEQ6M) and the dash-dotted (MSTW) lines and at the 
NLL by the difference between the continuous (CTEQ6M) and the dashed (MSTW) lines. The 
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Figure 5: The dependence of the total cross section pp gg on the renormalisation and fac- 
torisation scale, fx = fjLR = fip = ^nT-g, for r = rrig/mq = 1.2 and: (a) rrig = 0.5 TeV, (h) 
rrig = 1 TeV and (c) rrig = 2 TeV. The four lines in each plot correspond to the NLO and NLL 
cross section, each one evaluated with CTEQ6M and MSTW pdfs, see the legend. 
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Figure 6: The dependence of the total cross section pp qq on the renormalisation and fac- 
torisation scale, fi = fjLR = hf = (,n^q, for r = nig/mq = 1.2 and: (a) ruq = 0.5 TeV, (h) 
ruq = 1 TeV and (c) rriq = 2 TeV. The four lines in each plot correspond to the NLO and NLL 
cross section, each one evaluated with CTEQ6M and MSTW pdfs, see the legend. 
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absolute difference between tfie CTEQ6 curves and the MSTW curves is ratlier similar at tfie 
NLO and at tfie NLL accuracy. It is clearly visible in Fig. 5 that the uncertainty introduced by 
the pdfs grows with the increasing gluino mass. In particular, for nig = 0.5 TeV this uncertainty 
is smaller than 3%, for rUg = 1 TeV it is smaller than 5%, and for rUg = 2 TeV it reaches 
already about 25%. We rephrase that the probable reason for the strong dependence of the 
cross sections on the pdfs at large gluino mass is the fact that the currently available pdfs are 
poorly constrained at large scales and at large parton x. This uncertainty should be, however, 
substantially reduced after high pr jet measurements are perfomed at the LHC. As seen in Fig. 
5, the variation of the pp gg cross section with the scale ^ = ^rug is substantially reduced 
after inclusion of the soft gluon resummation over the whole ^ range. 

The case of squark-antisquark production is illustrated in Fig. 6. Clearly, for qq production, 
the soft gluon resummation introduces much smaller reduction in the scale dependence over the 
whole ^ range than it was in the case of the gg production. The relative uncertainty of the 
PP QQ cross section due to the choice of the pdfs varies from about 3% at rriq = 1 TeV by 
about 5% at niq = 0.5 TeV to about 8% at rriq = 2 TeV. Again, this uncertainty should be 
reduced after suitable measurements at the LHC are performed. 

We also study the corrections to the cross sections for pp gg and pp — > qq coming from the 
leading Coulomb exchanges beyond the NLO accuracy. We define the corresponding Coulomb 
if -factor 



IJ?. The obtained numerical results for the relative Coulomb correction, i^coui — 1, to the gg 
and qq cross sections are shown in Fig. 7a and Fig. 7b, respectively, and compared to the 
corresponding relative soft gluon correction, i^NLL — 1- For the gluino pair-production we find 
that the Coulomb if-factor exhibits only weak dependence on the gluino and squark masses. In 
more detail, i^coui — 1 takes values between 5% and 6%. In fact, the Coulomb corrections exceed 
the soft gluon corrections for rrig < 0.5 TeV and are relatively important up to rrig = 2 TeV. For 
the qq production -fCcoui ~ 1 takes the values between 2% and 3.5%, depending on rriq and rrig. 
The importance of the leading Coulomb corrections in qq production can be seen in Fig. 7b. 

The NLL soft gluon corrections and the Coulomb corrections beyond the NLO approxima- 
tion (57) can be combined additively. We define the if-factor that accounts for both types of 
corrections in the following way: 



where we use 



Eq. (57) to calculate 5a^ 



-K'Coul 



a 




-^NLL+Coul 




(p, m2,{/x2}) 



(60) 



20 




0.2 0.5 



m; 



1 1.5 
[TeV] 



a) 



0.2 0.5 



m; 



1 1.5 
[TeV] 



b) 



Figure 7: The relative corrections, i^NLL — 1 <ind -f^coui — 1; to the NLO cross sections for the gg 
(a) and the qq (b) production at the LHC as a function of gluino and squark mass, respectively; 
r = mg/ruq. 
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Figure 8: The K -factors, i^NLL+Coul for the gg (a) and the qq (b) total production cross sections 
at the LHC as a function of gluino and squark mass, respectively; r = nig/mq. 



21 



The obtained -fC-factors, i^^NLL+Coub in the gg and qq production at the LHC are shown in Fig. 
8a and Fig. 8b, respectively. The results shown in Fig. 8 are our most complete estimates of the 
higher order QCD corrections in these processes. 
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Figure 9: The K -factor, -ft^NLL' /^'^ partonic channels in the hadronic process pp gg: (a) 
qq gg and (h) gg ^ gg- 

It is interesting to investigate in more detail the soft gluon corrections in partonic subchannels 
of the gg and qq hadroproduction at the LHC. For this purpose we define the NLL i^T-factors 
restricted to the subprocesses occurring through the qq and the gg partonic collision, i.e. we 



define K^^^ 



(match) 



^pp^qq^kl/'^pp-Sqq^kl ^'^'^ "^NLL ~ '^pp^gg^kl 



(match) 



/<^pp^gg^kh where, as usual. 



kl = gg or kl = qq. The results for i^^LL ™ the gg and in the qq production are given in Fig. 9 and 
Fig. 10, respectively. The NLL i^-factors grow with increasing masses of the produced particles 
in all cases. This is expected since the importance of the threshold logarithms strengthens with 
the higher masses of the produced particles. 

It is clearly visible from Fig. 9 and Fig. 10 that the size of the soft gluon corrections increases 
with higher colour charges of all particles involved in the partonic reaction. The NLL correction 
for the gg mediated processes (see Fig. 9a and Fig. 10a) is much larger than the NLL correction 
for the qq mediated processes (Fig. 9b and Fig. 10b), -S^nll ~ 1 > -^nll ~ both for the gg and 
in qq production. More specifically, -f^^LL ~ 1 is 3 — 5 times larger than i^^LL ~ depending on 
the final state particles and their masses. By comparing the X-factors with the same partonic 
initial state in the production of gg and qq (i.e. by comparing Fig. 9a with Fig. 10a and Fig. 
9b with Fig. 10b) one finds that the soft gluon effects are larger for gg production than for qq 
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Figure 10: The K -factor, -R'nlL' /'^'^ partonic channels in the hadronic process pp qq: (a) 
qq qq and (h) gg qq. 



in both partonic subchannels. Such results are expected since due to larger colour factors the 
soft gluon corrections should be more pronounced when the particles with higher colour charges 
take part in the hard scattering. 

We also study the relative importance of the partonic subchannels in the considered sparticle 
production processes. For this purpose we define the fractions -^channel ~ ^pp^ij^ki / CTpp^ki for 
each total hadronic cross section, Upp^ku where kl = gg or kl = qq. The components app^ij^ki 
of the hadronic cross sections come from the contributions with intermediate partons ij. We 
shall denote: ij = qq {qq in the figures) for the combined qq and qq contributions, ij = gq for 
the combined gq, qg, gq and qg contributions, and ij = gg for the gg contribution. For the 
considered processes the ij = gq contribution does not appear at the LO. In Fig. 11 and Fig. 12 
we show the obtained values of -^channel 99 ^^^^ production at the LHC. The results 

are given at the LO (Fig. 11a and Fig. 12a), at the NLO (Fig. lib and Fig. 12b), and at the 
matched NLL accuracy (Fig. 11c and Fig. 12c). Clearly, the dominant contribution to the gg 
hadroproduction comes from the gg subprocess and the qq hadroproduction is dominated by the 
qq subprocesses except for the lower range of the squark masses where qq and gg contributions 
are similar. The inclusion of the NLO effects leads to a relative enhancement of the gg channel, 
and the values of -Rchannei ^'^ change visibly between the NLO and the NLL results. The 
gq channel appears beyond the LO and plays a minor role. Note that -^channel negative for 
the gluino-pair production (see Fig. lib and Fig. 11c), which is explained by the fact that this 
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Figure 11: Relative contributions of the partonic channels to the pp — > gg cross section at the LO 
(a), NLO (h) and NLL-matched accuracy. The thick lines correspond to r = 0.5, the medium 
lines to r = 1.2, and the thin ones to r = 2.0. The gg lines in (h) and (c) incidentally overlap 
for r = 0.5 and r = 2.0. 




Figure 12: Relative contributions of the partonic channels to the pp — > qq cross section at the 
LO (a), NLO (h) and NLL-matched (c) accuracy. The thick lines correspond to r = 0.5, the 
medium lines to r = 1.2, and the thin ones to r = 2.0. 



24 



contribution enters only beyond the LO where the simple probabilistic interpretation of the pdfs 
is lost. 

It is interesting to check what fraction of the NLO corrections is generated by the soft gluon 
contributions. In order to answer this question we truncate at the NLO the expansion of the 
resummed hadronic cross sections (Jp^^--^f^i mediated by the partonic subchannel ij and define: 



(exp) _ (res) 

pp^ij^kl pp^ij^kl 



. In Fig. 13 and in Fig. 14 we plot the obtained values of the 

(NLO) ^ Si' 



ratio, app]^-j_^^.i / a^^j^ki, for the gg production {kl = gg) and the qq {kl = qq) production, 
respectively. It follows from the figures that the soft gluon correction provides the dominant part 
(more than 2/3 in the studied cases) of the NLO correction in both the qq and gg subchannels, 
and both for the gg and qq production. The cusps visible in the curves in Fig. 13 originate from 
the NLO supersymmetric QCD correction to the pp — > gg cross section. In more detail, the 
stop-top loop contribution to the gluino self-energy exhibits the singular behaviour in vicinity 
of threshold for the decay of gluino into the top quark and its superpartner, the stop. 

Finally, wc analyse the effect of soft gluon resummation in the partonic colour channels. 
Recall that the soft gluon corrections at threshold do not lead to mixing of the s-channcl colour 
representations. Let us denote by cr'^p^^j^^i. j the contribution to the pp kl cross section at 
the LO coming from the partonic channels ij in the s-channel colour representation /, and by 
o'pp^ij^kl- 1 analogous contribution to the resummed hadronic cross section (not matched 
to the NLO). Furthermore, we define ratios 

^color - ^F^— (Ol) 

/ ,J ^pp—^ij^kl; I 

that correspond to the relative contribution of the colour representation / to the partonic channel 
ij of the hadronic cross section (we suppress here the indices of the initial, intermediate and 
final state particles in ^color)- -^^S- ^"^^ -^^S- show the ratios R^^^i^^ for gg and qq 

production, respectively, both for the LO cross sections and the resummed cross sections with 
r = 1.2. In general we find that the dependence of ^color °^ mass of produced particle is 
weak. Wc find that the colour singlet s-channel representation dominates in the qq channel of 
pp qq, and the colour octet s-channel representation dominates in the qq channel of pp gg 
and in the gg channel of pp qq. In all these channels the soft gluon effects only weakly affect 
-^color- ^^^'^ ^^^'^ of 99 channel of pp gg, the pattern is more interesting. Wc observe the 
following hierarchy of the contributions from the three colour channels, -Rcoior > -^coior > -^coior- 
The largest contribution coming from the 27 s-channel representation is most strongly enhanced 
due to soft gluon effects. 

In Fig. 17 and Fig. 18 we show the effect of soft gluon corrections on the cross sections in 
partonic-colour channels for the gg and qq production, respectively. Wc plot the i^-factors, K = 
af^^^/af''^^ = <^fj^-^kl I / "^ij^-^kl I these channels. The pattern is quite clear: the enhancement 
from soft gluon resummation for the gg partonic channels is strong and steeply rising with mass 
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Figure 13: The fraction of the NLO correction exhausted by the soft gluon contribution, 
c7(^''pVo-(nlo); for (a) qq^ gg and (b) gg gg; r = mg/rriq. 




Figure 14: The fraction of the NLO correction exhausted by the soft gluon contribution, 
cr(''''PVo"NLO, for (a) qq^ qq and (b) gg qq; r = rug/mg. 
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Figure 15: The decomposition of the cross section for the gg hadroproduction into the partonic 
colour channels, -Rcolor' /^'"■' ('^) ~^ 99 ^'^'^ (^) 99 ~^ 99' ~ l-^- thick lines correspond 
to the resummed cross sections and the thin ones to the LO cross sections. 
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Figure 16: The decomposition of the cross section for the qq hadroproduction into the partonic 
colour channels, -R^oior' f^''^- (^) W ~^ 99 ^^'^ (^) 99 ~^ t = The thick lines correspond 
to the NLL cross sections and the thin ones to the LO cross sections. 
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Figure 17: The K-factors, K = Gj 10^^*^^ for partonic colour channels: (a) qq ^ gg and (h) 
gg gg; r = 1.2. The dimensions of the s-channel colour representations are indicated in the 
legend. 
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Figure 18: The K-factors, K = aj^ /a^^^\ for partonic colour channels: qq qq and (h) 
gg — > qq; r = 1.2. The dimensions of the s-channel colour representations are indicated in the 
legend. 
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of the produced particles. This is not the case for the qq channel, where the soft gluons effects 
are small and weakly depending on the mass of the produced particles. Note that the if-factors 
in the qq channels are smaller than 1 for most values of the masses. In these channels, the 
size of the enhancement coming from soft gluon emission does not sufficiently compensate the 
suppresion due to changing from the LO to the NLO approximation of the pdfs and running 
CKg- For all considered channels we observe larger soft gluon enhancement for higher s-channel 
colour representations. 

7 Conclusions 

In this work we have studied the effect of soft gluon emission on the production of gluino-gluino 
and squark-antiquark pairs in the proton-proton collisions at the LHC. A detailed description 
of the derivation of the soft anomalous dimension matrices for the qq — ^ gg and gg — > gg 
scattering processes [13] has been given. It should be stressed that these matrices govern the 
soft non-collinear gluon radiation in the hadronic pair production of any heavy particles in the 
colour octet representation, like for instance the supersymmetric heavy colour scalar particles 
considered in [44] . In the threshold limit the obtained soft anomalous dimension matrices become 
diagonal in the s-channel colour basis. The diagonal elements of the soft matrices have been 
found to correspond to values of the quadratic Casimir operators of the SU(3) group for the 
outgoing pair of heavy particles. This observation confirms simple physical picture of the direct 
dependence of the soft gluon radiation at threshold on the total colour charge of the final state. 

The NLL resummation of the soft gluon corrections to the hadroproduction processes, pp 
qq and pp gg at the LHC has been explicitly performed for the values of masses of the produced 
particles between 0.2 TeV and 2 TeV and for the gluino-to-squark mass ratio 0.5 < r < 2. 
The obtained results have been matched to the corresponding cross sections computed at the 
NLO accTiracy. We have determined the NLO i^T- factors, -S'nll, for the squark-antisquark and 
gluino pair-production for these processes and studied their scale dependence. Futhermore, we 
have investigated the dependence of the NLL X-factors on the choice of the parton distribution 
functions. Additionally, the effect of the leading Coulomb corrections on the considered hadronic 
cross sections has been discussed. We have also analysed the partonic channel decomposition of 
the pp gg and pp qq total cross sections. In particular, we have studied the effect of the 
NLO and NLL corrections on the total cross sections in the partonic channels. Finally, we have 
decomposed the LO and NLL hadronic cross sections for the partonic subchannels in the colour 
basis and observed that the soft gluon enhancement grows with the the total colour charge of 
the pair of the produced particles. 
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A Mellin-moment transforms of the leading-order cross sections 

We define 

/■I , z^+^ , / 2r2(l + VT^) + (1 - r^)z\ 

Jn = / dz J- r^log' 



/o l+{^)z ''\2r\l-Vl^z) + {l-r^)zj ' 
Bn = /3(iV + 1,1/2) , 

G^^ = 2F1 I l,l/2,A^ + 5/2, 



GS^ = 2i^i(l,A^ + l,A^ + 5/2,-i!:^^j , 
where r = mg/rriq. With these definitions, the LO cross sections in iV-space read 



X 



81 (2Ar3 + 9iV2 + 13A^ + 6) (r^ + 1) 
-(r^ + l)A^(2Ar + 3) + 2{r^N + N + 2){N + 1)^^^ + (r^ + 1){N^ + 3N + 2)G^^ 



(62) 



kjHf 



+ 



gg^qq,l 



~(0) 
(J 



noq L9(4Ar2 + 16Ar + 15) 
2Sij ((r2 + l)(3r2 + N (2r2 - 1) - 2) - 2r\N + 1)G^^^) 
27 (2iV3 + 13 Ar2 + 27N + 18) (r^ - 1) 

as^nBN {N^ + 3N + A)nf 

ruq^ 96 (27V3 + 137V2 + 27iV + 18) ' 
ttsV^iv (5iV3 + 32iV2 + 717V + 



8 99^99,1 ' 



_2 



~(0) 

'"qq-*gg,l 



+ 32N^ + 71N + 68) nj 
96 (4iV4 + 36iV3 + 119iV2 + 171iV + 90) 



_ TT-BjV 1 



~{0) 

"'qq-*gg,8 



54(2Ar2 + 7iV + 6)(r2 + l) [(^^ + ^K-^iV - 12) 
+ (87V + 8)gS^^ + 2(r2 + l)(iV + 2)Gg) + (r^ + 1) (2Ar2 + 7iV + 6) ^ 



m-g 



27 (4iV4 + 3g^3 + ii9jv2 + i7iiv + 90) (r^ - 1) 
Mir"^ - DN'^ - (43 - 25r^)iV2 - (162 - 99r2)A^ - 



2(r2 + l)((r^ - l)iV'^ - (43 - 2br'^)N'^ - (162 - 99r2)A^ + 108r^ 
+ 4(2iV(r2 - 1) + 6r2 + 3) (2iV2 + TAT + 5) cj^^ 
- 7(r^ - 1) {2N^ + 15Ar2 + 37 AT + 30) 
+ fr^ - 1) f4Ar^ + 36Ar3 ^ ;L19Ar2 + 171 AT + 9o) ^ 



153) 



(63) 



(64) 



(65) 



(66) 
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as'^TTBN 9 {N^ + 9N^ + 20N + 14) 

(67) 



m-g'^ 64 {2m + 15iV3 + 4oiv2 + 45jv + 18) ' 
a^^TTBN 9 {N^ + IIN'^ + 30iV + 26) 



^^"^ rug^ 64 (4iV5 + 40iV4 + 155iV3 + 290iV2 + 261iV + 90) ' 

where we put = Ois- For the numerical evaluation of Jjv we separately use two forms of its 
expansion 

/ , 2 \ 2*;+l 

* = rT7^E TTTi E ,t\ 2 ^(^ + ^.* + "- + 3/2), (69) 

TO=0 \ / A;=0 ' 

and 



(68) 



, _ - (l-r^Y 1 - (-1)- 

^ l + r2^^j ^l + r2y 1 + m^ /3(A;+l,m-A;+l) ^ ^ 



B LL and NLL functions 

The expressions for the resummed factors, expanded up to NLL, are 

logAi(iV,4m^/x2) ^i^^fff) {boasifi^)logN) logiV + (to a,(M') log iV, 4771^, m') , (71) 
logAglfc,,,(iV,4m^/x2) (6oa.(/.^)log7v) (72) 

with 

5f^(A) = ^ [2A + (1 - 2A) log(l - 2A)] , (73) 



#(A,4m^M^) = -^log(l-2A) + ^ 



2A + log(l - 2A) + ^ log2(l - 2A) 



7r6o 27r6Q 

^ [2A + log(l-2A)l - ^,„g(l-2A)log(j;^) , (74) 



.(2) _ log(l^A) ^(1) 

where we took ix = jxp = jiR and 60 and hi are the first two coefficients of the QCD /^-function. 



, llCA-^TRUi . nCl-lQCATRnf-QCpTRni 
= 12^^ ' = 24;^^ • 

The values of the coefficients A^^ \ Af\ / defined in Section 3.1. 
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